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Abstract 

We present non-linear transport measurements on suspended, current annealed 
bilayer graphene devices. Using a multi-terminal geometry we demonstrate that 
devices tend to be inhomogeneous and host two different electronic phases next 
to each other. Both of these phases show gap- like features of different magnitude 
in non-linear transport at low charge carrier densities, as already observed in 
previous studies. Here, we investigate the magnetic field dependence and find 
that both features grow with increasing field, the smaller one with 0.6 meV/T, 
the larger one with a five to ten times higher field dependence. We attribute 
the larger of the two gaps to an interaction induced broken symmetry state and 
the smaller one to localization in the more disordered parts of the device. 

Keywords: A. Bilayer graphene; A. Multiterminal device; D. Broken 
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1. Introduction 

The isolation of monolayer graphene [l[ has started a new research field 
on two-dimensional carbon based electronic systems. Whereas the interest in 
monolayer graphene mainly stems from its unusual pseudo-relativistic electron 
dispersion, bilayer graphene promises to be interesting because of its strongly 
interacting charge carriers. Weak screening in suspended samples and small 
carrier densities at the charge neutrality point lead to large values for the inter- 
action parameter rg about one order of magnitude larger than in conventional 
2DEGs of comparable carrier density [||. For interaction effects to be observ- 
able in electronic transport experiments, bilayer graphene has to be very clean 
to achieve low charge carrier densities. To reach this regime, the normally used 
SiC>2 substrate has to be removed by underetching, leaving the bilayer freely 
suspended between the contacts [1, @| ■ The suspension allows to current anneal 
the bilayer graphene flake by Joule heating induced by passing a large current 
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through it. During this heating cycle adsorbates may (at least partially) evap- 
orate. 

The exact nature of the electronic groundstate of bilayer graphene has been 
investigated theoretically for several years d @, 0, H i, \m Ell, El El, El ■ The 



fourfold degeneracy that stems from the spin and valley degrees of freedom of 
the charge carriers will be broken by interactions, and new groundstates with 
partly or totally lifted degeneracies are expected to occur. Interaction induced 
broken-symmetry states that are being discussed as possible candidates for the 
groundstate of bilayer graphene include the quantum anomalous Hall insulator 
the quantum spin Hall insulator (QSH) 11211 , the nematic phase 
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the layer antiferromagnetic phase (LAF) [l_2j, |17|, [l_3( , and possibly more. 



Most, but not all of these states show a gap in the excitation spectrum around 
Ep=0, leading to a vanishing conductivity around the charge neutrality point 
(CNP). 

We can probe the electronic properties of the graphene bilayer as a func- 
tion of temperature, charge carrier density and magnetic field by measuring the 
differential conductance dl/dV — Gd as a function of a bias voltage V s d- We 
investigate low conducting and insulating samples with similar features as in a 
previous study, that exhibit two distinct gap- like features in Gd(V s <i)[lIj], and 
that mainly differ in the minimum conductance at V s d — 0. Here, we further 
investigate these kind of samples by looking at the magnetic field dependence 
of Gd{V s d) [lij] . Our results suggest that the samples are inhomogeneous after 
current annealing and are composed of two different electronic phases. A pos- 
sible reason for this is that the process of current annealing does not clean the 
samples in a uniform way, but creates almost disorder free areas in the center of 
the samples and leaves some disorder at the edges, which leads to a coexistence 
of electronic phases next to each other. 



2. Material and methods 

Bilayer graphene was deposited onto a p + -doped silicon wafer covered with 
300 nm silicon dioxide by micromechanical cleavage of natural graphite. After 
deposition, a thin (5 nm) film of aluminium was evaporated over the whole 
wafer. After exposure to air, this film oxidised and transformed to AI2O3 almost 
completely. Using standard electron-beam lithography, a PMMA etchmask was 
then structured for shaping the graphene flake using an Ar/02 plasma, and in a 
second lithography step, the electrode structures were defined. After each resist 
development, the exposed part of the oxidised aluminium film was removed in 
a 25% tetramethylammonium hydroxide (TMAH) solution, thereby eliminating 
any organic residues left on the graphene. The electrode structures (3 nm Cr and 
70 nm Au) were e-gun evaporated, and subsequently the graphene devices were 
annealed in high vacuum at mbar and 200 °C for several hours. The silicon 
dioxide was then partly removed in 2% buffered hydrofluoric acid, leaving the 
graphene suspended 160 nm above the substrate. The remaining silicon wafer 
acted as a gate with a voltage V g applied. The devices were mounted into a 
3 He cryostat in vacuum and subsequently current annealed at T = 1.5 K with 
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typically 0.4 mA per ftm of flake width [l8j]. For the present study, we fabricated 
devices with a four-terminal Hall cross geometry (insets of fig. [J); the longer 
axis was about 3 /xm long, the shorter 1.5 fxxn and the width of each individual 
arm was 1 jjm. Current annealing was done on pairs of contacts with the other 
two contacts floating. The annealing process was considered to be successful if 
a pronounced dip in G measured as a function of the gate-voltage V g appeared 
in the vicinity of V g = 0, i.e. at low charge carrier density of at least one contact 
pair combination. 

3. Results 
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Figure 1 : (a) Two-terminal conductance G of a Hall-cross device as function of the gate voltage 
V g at 1.5 K in device B2a. The possible configuration are shown in the inset. Measurements 
involving contact B show a pronounced dip in G at V g ~ -0.5 V. The remaining contact pairs 
have a flat gate response, (b) G as function of V g for all possible combinations of contacts 
in device B2b at 1.5 K. All curves show a minimum in G at V g ~ -0.5 V. The inset shows a 
scanning electron micrograph of the device. 

We present data taken on two devices that could successfully be current 
annealed. Fig. Q] summarises our findings. In the first device, B2a, the current 
annealing resulted in an inhomogeneous cleaning of the graphene. Fig. [T^,) 
demonstrates that only in two-terminal measurements of G(V g ) that involve 
contact B a CNP is visible with a minimum conductance of 0.9 e 2 /h. For all 
other contact pairs G only fluctuates around 3 e 2 /h as V g is changed and does 
not display a notable minimum. Further annealing rendered the sample more 
inhomogeneous, leading to gate characteristics with several CNPs at different 
V g . In contrast, in B2b the current annealing yielded a more homogeneous 
device. Fig. [Tb) shows G as function of V g for B2b. A CNP at V g w -0.5 V is 
well developed for all contact pairs. A minimum of G ~ 0.2 e 2 /h at the CNP 
can be found in G(V g ) measured between contacts A and B. At higher doping, 
the conductance saturates at 3-4 e 2 /h. 

Performing differential conductance measurements Gd{V s d) reveals further 
differences between B2a and B2b. In fig. a colour scale plot of Gd in B2a as 
function of V g and V s d at 230 mK is shown. The line profile at the CNP at V g w - 
0.5 V exhibits two gap-like features: a smaller one, called 5 in the following, is 
present up to \V s d\ = 0.6 mV and reduces Gd from 1 e 2 /h to 0.5 e 2 /h at zero 
V s d- The larger feature A is visible up to |V^<j| = 3.5 mV and lowers Gd from 
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Figure 2: (a) Differential conductance Gd as function of gate voltage V g and bias voltage 
V a d for device B2a at 230 mK. Curves overlaid in white are G £ j(V st j)-traces taken at gate 
voltages indicated by the black arrows. Close to the CNP at V g ss -0.5 V Gd is suppressed 
to below 1 e 2 /h for small bias voltages. At \V a d\ = 3.5 mV Gd reaches a maximum and then 
saturates around 3 e 2 /h. Away from the CNP, only a weak dependence on V s d is observable, 
(b) Temperature dependence of Gd as function of V a d of device B2a at the CNP. Two distinct 
gap features can be identified. The smaller one, 8/e fti 0.6 mV, is reduced when T is raised 
from 0.23 K to 1.5 K. The larger gap, A/e m 3.5 mV, remains unchanged, (c) Gd as function 
of V g and V sd for device B2b at 230 mK. At the CNP, V g S3 -0.5 V, and around zero V sd the 
device is insulating. Changing V g or applying V s d recovers Gd ~ 3 e 2 /h. (d) Temperature 
dependence of Gd as function of V s d at the CNP of B2b. The smaller gap, S/e m 1.5 mV, 
is reduced when T is increased from 0.23 K to 1.5 K. The larger gap, A/e PS 4 mV, is only 
reduced at 4 K. 

3 e 2 /h to about 1 e 2 /h at zero V s d- Moreover, A shows a BCS-like overshoot 
in Gd- Moving V g away from the CNP rapidly suppresses A. On the other 
hand, 8 persists as a small zero bias anomaly even in the metallic regime and is 
visible in the whole accessible gate voltage range. In fig. [2J)) the temperature 
dependence of Gd at the CNP of B2a is plotted. The smaller gap feature, S, 
is totally suppressed if T increases from 0.23 K to 1.5 K. The larger feature 
A remains unchanged. The temperature dependence G(T) measured at large 
carrier density can be fit with a log(T) dependence for T<1K. 

Device B2b differs from B2a by having a full gap, which reduces G to zero 
at the CNP. Fig. [2b) shows G d as function of V g and V sd at 230 mK. The 
line profile at the CNP is qualitatively similar to B2a, but Gd is completely 
suppressed for small V s d- Furthermore, the small gap feature 5 is not readily 
visible but manifests only as a change in the slope of G d at \V sd \ =1.5 mV. 
The larger feature A is shifted to higher voltages as well, with \V 3 d\ = 4 mV, 
and shows a larger overshoot before Gd saturates at 2.3 e 2 /h. Tuning V g away 
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from the CNP quenches the larger feature A but leaves the smaller 6 intact in 
agreement to what is observed in B2a. In fig. [2jl) the temperature dependence 
makes 5 better visible. The small gap feature shrinks as T is raised from 0.23 K 
to 1.5 K, but still reduces Gd to zero at V s d — mV. A further increase to 
4 K completely suppresses <5. The temperature change from 0.23 K to 1.5 K 
only reduces the overshoot in the larger feature A. Increasing T to 4 K yields 
a reduction of A and removes the overshoot. 
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Figure 3: (a) Two-terminal conductance G as function of bias voltage V a d and perpendicular 
magnetic field B of device B2a at the CNP and 230 mK. The device is low-conducting at 
B = T and smoothly evolves into an insulating state for B > 0.5 T. (b) The corresponding 
line profiles. The dashed lines are guides to the eye for the linear scaling of A and S with B. 
A scales as 2.8 meV/T, whereas <5 grows with 0.6 meV/T. (c) G as function of V s d and B of 
device B2b at the CNP and 230 mK. The device is insulating even at B = T. The gap A 
increases linearly with B. (d) The corresponding line profiles. The dashed line is a guide to 
the eye for the evolution of A with 6 meV/T. 

To find out about the origin of the two different features S and A we perform 
differential conductance measurements Gd{V s d) as a function of perpendicular 
magnetic field, shown in figure |31 Both S and A increase linearly with \B\ 
but with different slope. Taking the two inflection points of the Gd{V s d) traces 
as a measure for S and A, we can determine 85/ dB to 0.6 meV/T, which is 
considerably larger than the Zeeman energy of a spin 1/2 (0.12 meV/T for 
g=2). Looking at the large gap A, we find dA/dB — 2.8 meV/T for sample 
B2a (Fig. |3b) and 6 meV/T for sample B2b (Fig EH). As sample B2b showed 
a full gap already at B = 0, the magnetic field dependence of S could not be 
determined here. Similarly, Velasco et al. observed 5.5 meV/T 19] for a fully 
gapped bilaycr graphene device for a gap of similar size as A. 
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The experiments clearly show two distinct gap-like features that occur in 
different ranges of charge carrier density and temperature, and that show a 
clearly different magnetic field dependence. We think that this observation 
can only be explained by the coexistence of different electronic phases next 
to each other in different parts of the samples. This conjecture is supported 
by the results of the current annealing procedures shown in figure [TJ that can 
be understood by the existence of almost disorder-free, clean regions that are 
surrounded by more disordered areas where the current annealing has been 
less efficient. In the disordered regions the CNP lies outside the charge carrier 
densities that are accessible to us, due to chemical doping by adsor bates and 
residues. In this picture, sample B2a would be composed of a clean region close 
to contact B, and a more disordered rest of the sample. In sample B2b the clean 
region seems to be larger, occupying the center of the sample, but different CNP 
conductances between different terminals indicate that even in sample B2b there 
are regions with stronger disorder that are not evenly distributed, as illustrated 
in Fig. gj 

We assume that the clean phase is showing an interaction induced gap A 
close to the CNP with vanishing conductivity, and a distinct, BCS like shape. 
The large increase of A with B of several meV per tesla has been predicted by 
theory for certain broken symmetry states like the QAH or QSH [l2|, Ell ■ 

a) B = 0T b) B>0T 




Figure 4: Sketch of the graphene device with the clean (light grey) and the disordered regions 
(dark grey) at the charge neutrality point, (a) In the absence of an external magnetic field, 
the clean region is gapped with A and the contacts do not couple to it. The disordered region 
remains conductive (red arrows), but experiences a transport gap 6 due to localization, (b) In 
the presence of a perpendicular magnetic field the clean region expands, leading to a narrowing 
of the conducting areas <5 and to a suppression of conductance. 



The total conductance of sample B2a is composed of two parallel conduction 
channels, one through the clean region with gap A in the center of the sample, 
and another, more disordered one along the sample edges with gap- like feature S. 
The origin of the smaller feature S is probably disorder induced localization (20I 



2lLl22j. For a disordered device such as a graphene ribbon with a predetermined 



geometry the application of a magnetic flux through the sample should lift the 
localization at least in part. In contrast to this scenario, the small gap in 
our experiment even increases when a perpendicular magnetic field is applied. 
We explain this by a change in geometry. The application of a perpendicular 
magnetic field increases the size of the inner clean region and thereby narrows 
down the conduction channel at the sample edge (Fig. |4p) , leading to increased 
localization in the disordered part. The different response of the two conduction 
channels to magnetic field is reflected in the different magnetic field dependence 
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of dS/dB =0.6 meV/T and dA/dB=2.8 meV/T. 

In sample B2b the initial distribution of the two phases is such that already 
at B=0 the sample is mostly in the clean phase, and there is only a negligible 
conductance contribution of the disordered phase at the CNP. This is also the 
reason for the fact that only a large gap A can be observed in Gd(V s d), and that 
5 is feeble. 

Because the metallic contacts will act as a heat sink during current annealing, 
we also have to assume that the disordered phase persists close to them. In fact, 
a narrow area of organic deposits can be clearly identified close to the metal 
contacts in scanning electron micrographs. This disordered region will then 
decouple the clean area in the sample center from the contacts and thus prevent 
edge states, that exist in the case of a QAH or QSH phase to couple to the 
metal contacts. Consequently, any phase in the center of the sample would show 
insulating behaviour in two terminal conductance, no matter if it has edge-states 
(QAH, QSH) or not (LAF). As we have shown the current annealing process 
to be inhomogeneous, we think that 5 is a feature of the disordered graphene 
regions that persist at the edges of the device. However, as the size of the 
disordered region will vary from sample to sample, any distribution of 8 versus 
A -gaps may be realized. Different distributions of clean and disordered regions 
in the samples will lead to different Gd(V s d) traces around the CNP, with a full 
suppression of conductance at V s d—0 when the clean region dominates, and some 
small, finite conductance when there is a considerable region of the disordered 
phase remaining. The application of a perpendicular magnetic field will always 
lead to a growth of the gapped A-phase at the expense of the conducting 8- 
phase, in this way opening a full gap in all samples at some finite value of B. 
We note that, although not discussed in detail, a small gap 8 was also observed 
in [l9| . It persists at high charge carrier densities and shows a small magnetic 
field dependence, in contrast to the stronger magnetic field scaling of A and in 
full agreement to our observations. 

4. Conclusions 

In conclusion, we measured nonlinear conductance on multi-terminal, sus- 
pended bilayer graphene devices in a Hall-cross geometry as a function of charge 
carrier density and perpendicular magnetic field. We find that current annealing 
cleans the graphene inhomogeneously, which leads to samples that are composed 
of two different regions with distinct electronic properties. Both of these regions 
show gap- like features at the CNP, as visible in differential conductance Gd{V s d)- 
We extract two gaps of different magnitude, that show a different temperature 
dependence and respond differently to magnetic field. Attributing the two gaps 
to clean and disordered sample regions, we find that the intial distribution of 
the two regions differs from sample to sample, but that the overall density and 
magnetic field dependence of the conductance of these two regions is similar. 
We argue that one of the two regions is a conductive, disordered phase local- 
ized at the sample edges that shows a zero bias anomaly over a large range in 
charge carrier density due to disorder induced localization, and a clean phase 
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in the center of the sample that shows an interaction induced gap close to the 
charge neutrality point with strongly suppressed conductance. The application 
of a perpendicular magnetic field favors the expansion of the clean phase in 
the sample center and shrinks the disordered conducting areas, and therefore 
pronounces the signature of the interaction effects in differential conductance. 
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